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The main aim of this research article is a parametric demonstration of irreversible Stirling cryogenic 
refrigerator cycles that includes irreversibilities such as external and internal irreversibilities. In addition, 
through this study, finite heat capacities of external reservoirs are considered accordingly. To reach the 
addressed goal of this research, three objective functions that include the input power of the Stirling 
refrigerator, the coefficient of performance (COP) and cooling load (R L ) have been involved in optimiza¬ 
tion process simultaneously. The first aforementioned objective function has to minimize; the rest objec¬ 
tive functions, on the other hand, have to maximize in parallel optimization process. Developed multi 
objective evolutionary approaches (MOEAs) based on NSGA-II algorithm is implemented throughout this 
work. Moreover, cold-side’s effectiveness of the heat exchanger, hot-side’s effectiveness of the heat 
exchanger, heat source’s heat capacitance rate, heat sink’s capacitance rate, temperature ratio , tem¬ 
perature of cold side are assigned as decision variables for decision making procedure. To gain a robust 
decision, different decision making approaches that include TOPSIS, LINMAP and fuzzy Bellman-Zadeh 
are used. Pareto optimal frontier was determined precisely and then three final outputs have been gained 
by means of the mentioned decision making approaches. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Stirling engine producing low temperature stream, is one of the 
crucial refrigeration cycle approaches. More description about 
coefficient of performance (COP) of the refrigeration systems are 
presented in Ref. [1]. Stirling engine producing low temperature 
stream also have relevant outputs [1] for the COP of the refrigera¬ 
tion systems. 

Leff and Teeters [2] have mentioned that since there is no ‘nor¬ 
mal peak’ in these systems, the straight-forward Curzon-Ahlborn 
[1] determination will not cover for the reversed Carnot cycle. 
The Lagrangian approach of unknown multiplier has been utilized 
by Blanchard [3] to determine the endoreversible Carnot heat 
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pump’s COP operating at lowest input power for a specific heating 
load. Recently, the efficiency of other refrigeration cycles such as 
Carnot, Brayton, Stirling and Ericsson has been investigated, by 
employing the context of the finite time thermodynamics 
approach [4-9 , the criteria of ecological [10-17] and thermo- 
economic model [18-21 for various working situations. 

Finite time thermodynamic (FTT) is one of the basic methods to 
analyze thermodynamic systems, though the Carnot cycle’s endo¬ 
reversible approach has been presented and progressed for 
20 years from 1950s and 1970s [22 . The FTT approaches are 
utilized through the actual apparatus and processes optimization, 
related to finite-time and finite-size limitations [23 . In recent 
years, utilizing FTT concepts for the optimization and evaluate of 
heat engines and refrigeration systems for various objectives opti¬ 
mizations have made considerable advance [24-30 . Hints for the 
optimum design of different actual refrigerators have prepared 
by these developments. 
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Nomenclature 


A area 

b\ related to Eqs. (15)-( 17) 

C heat capacitance 

C v molar specific heat capacity, J mol -1 K” 1 

d Index of deviation 

d i+ deviation of ith solution from the ideal solution 

d,_ deviation of ith solution from the non-ideal solution 

COP coefficient of performance 

Fi related to Eqs. (15)—(17) 

i ith objective 

j jth solution 

n working fluid’s mole number, mol 

p power, W 

Q heat transfer, J 

R universal gas constant, J mol' 1 K _1 

T temperature, K 

t cyclic period, s 

U overall heat transfer coefficient 


x temperature ratio of the Stirling cryogenic refrigerator 

X decision variable’s vector 

Greek letter 

A ratio of volume through the regenerative processes 

(volumetric ratio) 

s H hot-side heat exchanger’s effectiveness 

s R regenerator’s effectiveness 

s L cold-side heat exchanger’s effectiveness 

Subscripts 

c sink side/cold side 

H, h heat source/hot side 

L cooling/heat sink 

R regenerator 

s isentropic/ideal 

I, 2,3,4 states points 


To involve all of the above-proposed objects three objective 
functions counting input power, the coefficient of performance 
(COP) and cooling load (R L ) have been investigated. Furthermore, 
Six decision parameters that contains the heat capacitance rates 
of heat sink and heat source, the cold- and hot-side heat exchang¬ 
ers effectiveness, ratio of temperatures 0 a) and temperature of 
sink side (T c ) have been used in the multi-objective optimization 
process. 

Multi-objective optimization has extreme presence in different 
engineering problems such as skyline computation, vehicle routing 
issues and so forth [31 -33 . Simultaneous consent of various objec¬ 
tive functions through for gaining optimum solution through the 
multi-objective optimization is a highly difficult target. The back¬ 
ground of the Evolutionary algorithms (EA) referred to 18th cen¬ 
tury in a stochastically effort to gain a solution of this type while 
this approach was evolved and applied gradually [34 . An accept¬ 
able route to a multi-objective quandary is summarized to prepare 
an assortment of routes, each of them consents the objectives at a 
reasonable degree without being overcome by the rest solutions 
[35]. Generally, Multi-objective optimization issues represent an 
uncounted collection of outputs which is called Pareto frontier, 
who examined vectors illustrating the prior feasible solution in 
the space of the objective function. Multi-objective optimization 
of system of energies and various thermodynamic have been the 
focus of many studies these days [36-43 . Recently, Ahmadi et al. 
[44,45] have developed a multi-objective optimization method 
for designing a power Stirling engine by applying the Genetic Algo¬ 
rithm (GA) which is known as a branch of evolutionary algorithm. 

In the current study, multi-objective optimization algorithms 
have been executed to maximize the coefficient of performance 
(COP) and cooling load (R L ) and minimize the input power. To ob¬ 
tain the optimum solution, three robust decision making methods 
were used and error analysis of the outcomes gained from each 
decision making techniques was employed to determine effective¬ 
ness and accuracy of each optimum route that obtained from the 
mentioned decision making methods. 

2. System description 


and two isothermal processes as displayed schematically in 
Fig. 1. The expansion stroke of an actual cycle can be determined 
in the addressed cycle as an isothermal process 1-2 join irrevers¬ 
ible isothermal heat addition owing to the finite heat capacity’s 
heat source that temperature differ from T u to T L2 at temperature 
T c . Adding of the heat into the working fluid from the regenerator is 
assumed as an isochoric process like as process 2-3. Isothermal 
processes 3-4 assumed for the compression stroke join irreversible 
heat declining to the finite heat capacity’s heat sink that tempera¬ 
ture changes in interval of T m to T H2 at temperature T h . Moreover, 
isochoric process 4-1 represents the heat declination from the 
working fluid to the regenerator. As mentioned earlier, the heat 
transfer in the mentioned processes such as 1-2 and 3-4 in an 
actual cycle have to happen in finite time. To reach this, these 
mentioned heat processes have to be preceded throughout a finite 
temperature variation and consequently introduced as external 
irreversibilities. 


3. Thermodynamic analysis 

The value of the absorbed heat of the cycle at temperature T c 
from the source represent by Qc and the value of the released heat 
of the cycle at temperature T h from the sink represent Q/ 7 therefore: 



Magnetic material or a gas can play role of working fluids in the 
Stirling cycle while various working fluids have various effective¬ 
ness properties. If the implemented working fluids assumed ideal 
gas in the cycle the referred cycle separates into two isochoric 


S (kJ kg 1 K 1 ) 




Fig. 1 . Qualitative scheme of an irreversible Stirling refrigerator cycle [9]. 
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Qh = T h AS = C h (T h 2 - T m )t H = U H A H (LMTD) H t H (1) 

Q.c = T C AS = C L (T L1 - T L2 )t L = U L A L (lMTD) L t L (2) 

where entropy change denotes AS = nR In X, and X stands for the 
volume ratio and n denotes the mole number of the working fluid. 
C L ; C H and t L \ t H refer to the source/sink reservoirs heat capacitance 
rates and times of heat addition/rejection, HI and LI represent hot- 
side and cold-side, correspondingly. U H A H and U L A L indicate the 
products of area and overall heat transfer coefficient and ( LMTD) H 
and ( LMTD) l denote the log mean of the difference between temper¬ 
atures on sink-side and source-side, correspondingly, expressed as: 


( LMTD) h = 


(Th - Th i) - (T h - T H2 ) 


In 


( (Th~Tm) \ 

\T h -T H2 )J 



Or 


Qc i Qh /1 a \ 

= ( 14 ) 

h 1 h 

In which <p denotes the irreversibility factor and is lower than 
one for an actual cycle. So, from Eqs. (6)-(14), following expres¬ 
sions can be developed: 


X- (j) 


ShC h (xT c -T H i ; 


+ (slCUTu-Tc)) +M* V) 



R, 


(p-F i(x-l) 


£ H C H (xTc-T m ) 


+ 




£ ; lQTli -T c 


+b,(x- 1) 



(LMTD) l 


(Tli - Tc) - (T L2 - T c ) 


In 


ATii-TcA 

{(Tu-Tc)J 



Generally, these cycles do not states the situation of complete 
regeneration, So it is favorable to involve the losses of regenerator 
that is related to the difference of temperatures between the two 
mentioned isothermal processes as supposed by previous re¬ 
searches in Refs. [14-16], 

AQ R = nC v (\-s R )(T h -T c ) (5) 

n denotes the mole of the working fluid. C v denotes the working 
fluid’s specific heat capacity through the processes of the regenera¬ 
tive in unit of mole, s R represents the regenerator effectiveness’s. 

Thanks to the irreversibility impact of the finite heat transfer, 
the regenerative time have to finite, alike to the above mentioned 
isothermal processes [12-16], and comes to following expression: 


t R = t 3 + t 4 = 2<x(T h - T c ) 



In which a stands for the proportionality constant, is not de¬ 
pend on the difference of temperatures; however, relates on the 
characteristics of the regenerative mater. Therefore, the total cycle 
time t C y C i e expressed as follow [12-16]: 

tcycle — (£h + h + t R ) (7) 

To take into account all of the above mentioned irreversibilities, 
the net value of the released heat to the sink and absorbed heat 
from the source are determined by following expressions [9,48]: 

Qh = Qh ~ AQ r (8) 

Ql = Qc- AQ r (9) 

As known for any experts of engine systems, the power input, 
cooling load and the COP are the most crucial variables of any 
refrigerator. By executing the above formulations, following for¬ 
mulas can be developed for the power input, cooling load, and 
cooling COP as following as [9,48]: 


p _(Qh- Ql) 

^cycle 

(10) 

r l = Ql 

l-cycle 

(ii) 

C0P L = | 

(12) 


Dedicated to the thermodynamics second law for irreversible 
cycles, following formulation should be considered 


£-£<0 (13) 


8 

-e- 

X *Ti 

1 l "" k 

-e- * 

(17) 

T h 

X = T 

1 c 

(18) 

, 2a 

(19) 

1 nR In X 

P C„( 1-£r) 

Fl = R In A 

(20) 


4 . Multi-objective optimization with EA 

4A. Optimization via EA 

Genetic algorithm was developed by Holland (1960) in order to 
involve the mechanisms of biological evolvement into computer 
code and further optimization process [37,46 . Genetic algorithms 
(GAs) handle a repetitive, stochastic search approach to determine 
an optimal route and mimic in a simplified manner principle of 
biological evolution (Fig. 2). In this study, GA is subset of evolution¬ 
ary approach is used to obtain the Pareto frontier. This approach is 
a robust optimization method for nonlinear issues [34-37,46]. To 
determine a ranking between the obtained routes of a multi-objec¬ 
tive optimization issue, optimization of Pareto is the vital key, till 
to figure out whether a route is actually one of the prior feasible 
solutions [46]. NSGA-II algorithm as a multi objective optimizer 
has been executed till to figure out the Pareto frontier from the 
GA. In Refs. [34-37,40] more about evolutionary algorithms have 
been presented. 

4.2. Objective functions, limitations and decision parameters 

The main objectives of the current communication are the sys¬ 
tem’s input power (P), the coefficient of performance (COP), and 
the cooling load (P L ) denoted by Eqs. (15)—(17), correspondingly. 

Following decision parameters have been involved in the 
current communication: 

T c : Temperature of cold side 
x: Temperature ratio ^ 

C H : Heat source’s capacitance rate 
C L : Heat sink’s capacitance rate 
s L : Effectiveness of the cold-side heat exchanger 
s H : Effectiveness of the hot-side heat exchanger 

However the above suggested decision parameters can be chan¬ 
ged throughout the optimization process, it should be varied with¬ 
in the practical or feasible range. Owing to this fact, following 
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Fig. 2. Chart Box of the multi-objective optimization executed in the current communication [36,37,40]. 


limitations have been considered for solving the referred objective 
functions [8,9,13,16,48]: 

0.5 Ch ^ 0.8 

(21) 

0.5 < s L ^ 0.8 

(22) 

235 ^ T c ^ 240 K 

(23) 

600 < C H ^ 1000 

(24) 

600 < C L < 1000 

(25) 

1.3 ^ x ^ 1.4 

(26) 


F- = —— fti , to minimize objectives (28b) 

lJ max(Fy) - min(Fy) v y 

In the current study, robust, familiar and usual type of decision¬ 
making techniques that includes TOPSIS, LINMAP and fuzzy 
Bellman-Zadeh, techniques are applied to gain final optimum 
solution is selected using these three methods. The LINMAP and 
TOPSIS methods utilize Euclidean non-dimensionalization while 
the Bellman-Zadeh technique executes the fuzzy non-dimension- 
alization and the following sections explain the mentioned 
decision-making techniques. 


5. Decision-making through the multi-objective optimization 

Process of making a decision for indicating the prior optimum 
route from practical outputs is needed through multi-objective 
optimization process. There are a number of methods for process 
of making decision in decision issues. To determine a prior opti¬ 
mum output from the Pareto frontier, these methods can be uti¬ 
lized. Before making any decision, objective space scales and 
dimensions have to unify while dimensions of different objectives 
might be different through multi-objective optimization issue. 
Accordingly, objectives vectors would be non-dimensioned before 
the decision-making process. Various techniques of non-dimen¬ 
sionalization employed for making a decision such as Euclidean 
non-dimensionalization, and fuzzy non-dimensionalization. 

5.2. Non-dimensionalization methods 

5.2.2. Euclidean non-dimensionalization 

If the objectives matrix for different solutions of the Pareto 
frontier is represented by Fjj where i is stands for the index for each 
route on the Pareto frontier and j denotes the index for mentioned 
objectives in objective space. In this method, a non-dimensional- 
ized objective, F”, is represented as following as: 

F 

Fjj = lJ .. . = for minimizing and maximizing objectives 

\/ZM 2 


5.2. Decision making techniques 

5.2.1. Bellman-Zadeh decision-making technique 

Prior optimum solution defined by Bellman-Zadeh approach 
can be indicated by intersection point between restrictions and 
fuzzy criteria meantime suggested by corresponding membership 
function. In Bellman-Zadeh technique, a membership function’s 
matrix is determined in which column illustrates the objective 
fuzzy membership function. In addition, rows of the mentioned 
matrix present values of the membership functions for each ob¬ 
tained outcomes from the Pareto frontier. Put it in other way, the 
number of matrix columns represents the number of proposed 
objectives and the number of matrix rows represents the number 
of obtained outcomes relied on the Pareto optimum frontier. Detail 
for method of definition for the membership function has been sta¬ 
ted in Refs. [36,37,47 . For the proposed solution, the fuzzy mem¬ 
bership function for each outcome is a group at the lower value 
limit of the all objectives membership functions. In this regard, a 
membership vector is determined who its components are the 
least objectives membership function for each outcome. Totally, 
the highest value of these lower limits i.e. a component of the 
membership vector with highest magnitude of membership func¬ 
tion is considered as a great optimum outcome. Let me rephrase 
that, throughout making a decision by executing Bellman-Zadeh 
fuzzy technique, an upper limit of lowest bound of objectives 
membership function is considered as the prior optimum solution. 
In Refs. [36,37,47] more details are presented. 


5.2.2. Fuzzy non-dimensionalization 

Throughout this technique, a non-dimensioned objective, Fjj, is 
expressed as following: 



— ^ ij ™ n (ft/) t0 maximize objectives 
max(Fy) - min(Fy) 


(28a) 


5.2.2. LINMAP decision-making technique 

A criterion of LINMAP decision making technique should be 
defined here that called “ideal point” who addressed to the data 
point on the Pareto frontier which each mentioned objective is 
optimized either singly. It would be stressed here that it is impos¬ 
sible to gain optimum objectives by utilizing single objective opti¬ 
mization because we have different objectives that optimized 
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simultaneously. Due to this fact, location of the ideal point is not on 
the Pareto frontier. LINMAP technique selects the optimum 
solution based on the distances of output on the Pareto frontier 
from the ideal point. In other words, the solution that has lowest 
deviation from the ideal point is considered as an optimum route 
through executing LINMAP method. More explanations of the 
LINMAP technique can be found in previous works of Ahmadi 
and colleagues [36,37,47]. 

5.2.3. TOPSIS decision-making technique 

Throughout TOPSIS technique another criterion should be 
added which defined as “non-ideal point” that refers to worst val¬ 
ues of the proposed objectives in the objective area. Keep in mind 
that same as previous decision making methods, ideal point is 
another important criterion for selecting optimum solution. So, 
distances between the solutions gained from the ideal and 
non-ideal data points are the selection criteria for determining 
optimum solution through executing TOPSIS method. Owing to 
this fact, throughout the TOPSIS technique the optimum solution 
has a maximum distance from non-ideal point and on the hand 
has a minimum deviation from ideal point. Further descriptions 
about TOPSIS technique have been given in Ref. [47]. 


6. Result and discussion 

The coefficient of performance (COP) and the cooling load (P L ), 
are maximized simultaneously and the input power is simulta¬ 
neously minimized by executing optimization approach that 
performs according to the NSGA-II procedure. 

In this respect, multi objective optimization is executing with 
the previously mentioned objective functions whose are formu¬ 
lated in Eqs. (15)-(17) limitations that are proposed with Eqs. 
(21)-(26). 

Decision parameters through the optimization algorithm are the 
cold side temperature, temperature ratio (j^), heat source capaci¬ 
tance rate, heat sink capacitance rate, effectiveness of the cold-side 
heat exchanger, effectiveness of the hot-side heat exchanger. 

Following Stirling cryogenic refrigerator specifications are 
considered to have same standard with previous literature [9], 
T„i = 300 K, T n = 250 K, s R = 0.9, 2 = 2, n = 1, R = 4.3(J mol” 1 K -1 ), 
C v = 15(J mol” 1 K _1 ), 0 = 0.8, a = 10' 5 (s K" 1 ). 

Pareto optimal frontier for three mentioned objective functions 
that include input power, cooling load and the coefficient of perfor¬ 
mance are represented in Fig. 3. Besides, selected data points 
thanks to the decision making approaches are depicted in Fig. 3. 



Fig. 3. Pareto optimal frontier in objectives’ space. 



Fig. 4. Pareto frontier (Pareto optimal solutions) for the input power versus cooling load using NSGA-II. 
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As it is shown Fig. 3 optimum results of TOPSIS and LINAMP 
methods are very close together and are more suitable than 
FUZZY method results, due to less power input and higher COP. 
Power input and COP are more effective than RL, on system 
performance. 

Pareto optimum frontier for cooling load and input power 
which are two objective functions are shown in Fig. 4. In addition, 
selected data points thanks to the decision making approaches are 
depicted in Fig. 4. As the system is more sensitive to the amount of 
input power and this sensibility is less to RL, results of TOPSIS and 
LINAMP methods are better than FUZZY method results. One of the 
optimization goals is to minimize the power input, which it leads 
to better performance of the system. 

Pareto optimum frontier for input power and the coefficient of 
performance whose are two objective functions are shown in Fig. 5. 
Moreover, selected data points thanks to the decision making ap¬ 
proaches are depicted in Fig. 5. The less input power and the higher 
COP lead to better system performance. According to this, results of 


TOPSIS and LINAMP methods are optimum in comparison with 
FUZZY method results which is displayed in Fig. 5. 

Pareto optimum frontier for cooling load and the coefficient of 
performance that are two objective functions are shown in Fig. 6. 
Furthermore, selected data points thanks to the decision making 
approaches are depicted in Fig. 6. 

Figs. 7-12 demonstrate the scattered distribution of the design 
variables. These plots can provide better vision of variation of 
the design variables from the Pareto frontier. The distribution of 
C H for the optimum data on Pareto front is demonstrated in 
Fig. 7 meanwhile the minimum and maximum boundaries of the 
variables are illustrated. Furthermore, Fig. 8 represents the vari¬ 
ables distribution for the optimum data on Pareto front for C L . In 
addition, the upper and lower boundaries of the addressed param¬ 
eters are illustrated. Also, the variation of thee addressed parame¬ 
ters for the optimal points on Pareto front for s H over population 
are depicted in Fig. 9 whereas upper and lower limit of s H referred 
to 0.5 and 0.8, respectively. The same analogy was employed for e Lf 



Fig. 5. Pareto frontier (Pareto optimal solutions) for the input power versus the coefficient of performance using NSGA-II. 
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Fig. 6. Pareto frontier (Pareto optimal solutions) for the cooling load versus the coefficient of performance using NSGA-II. 
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Population 

Fig. 7. The variation of C H on the Pareto front for the optimum points. 
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Fig. 8. The variation of C L on Pareto front for the optimum points. 
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Fig. 9. The variation of e h on Pareto front for the optimum points 
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Fig. 10. The variation of s L on Pareto front for the optimum points. 
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Fig. 11. The variation of T c on Pareto front for the optimum points. 
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Fig. 12. The variation of x on Pareto front for the optimum points. 
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Table 1 

Comparison of the final optimum solutions of this work (multi objective-multi variable) with relevant results gained in Refs. 1 9,48] (single objective-single variable optimization) 
(■ n=l , C v = 15 (J mol -1 K _1 ), R = 4.3 (J mol” 1 K _1 ), X = 2). 


Decision making method 

Decision variables 





Objectives 



£h 


T c 00 

X 

C H (W/K) 

Cl (W/K) 

P(W) 

Rl (W) 

COP 

TOPSIS 

0.79427 

0.79931 

239.784 

1.30002 

986.888 

990.969 

2067.701 

2835.817 

1.371483 

LINMAP 

0.79671 

0.79758 

239.849 

1.30116 

987.897 

990.986 

2105.365 

2879.062 

1.367488 

Fuzzy 

0.79919 

0.79919 

239.737 

1.30811 

987.796 

990.983 

2290.827 

3077.963 

1.343604 

Ref [9] 

0.8 

0.8 

240.5 

1.40 

1000.0 

1000.0 

4423.765 

4414.133 

0.998 

Ref [48] 

0.8 

0.8 

- 

1.90 

1000.0 

1000.0 

2710 

2370 

0.875 


Table 2 










Error analysis of the results gained from the three aforementioned decision making approaches in the basis of MAPE. 





Decision making .method 

TOPSIS 



LINMAP 



Fuzzy 



Objectives 

P 

Rl 

COP 

P 

Rl 

COP 

P 

Rl 

COP 

Max error % 

6.77 

11.78 

5.51 

5.43 

10.47 

5.81 

5.66 

10.35 

5.56 

Average error % 

2.78 

8.07 

5.44 

2.70 

7.62 

5.41 

2.81 

8.01 

5.36 


dedicated to this fact, the assortment of the mentioned parameters 
for the optimum data on Pareto front for the addressed parameter 
over population is illustrated in Fig. 10 and upper and lower 
restrictions of s L are 0.5 and 0.8, correspondingly. 

Fig. 11 depicts the variation of T c for the optimum points on Par¬ 
eto front though the lower and upper magnitudes of the referred 
parameters are marked. Finally, variables distributions for the opti¬ 
mum points on Pareto front for x over population are exhibited 
through Fig. 12 and upper and lower ranges of x eluded to 1.3 
and 1.4, respectively. 

The optimal outputs obtained from three mentioned decision 
making methods (Bellman-Zadeh, TOPSIS, and LINMAP) that in¬ 
clude decision parameters and objective functions are demon¬ 
strated throughout Table 1. 

Table 1 compares solutions gained from our suggested decision 
making techniques with relevant solutions found in Ref. [9,48]. 
Where Table 1 demonstrated outputs gained from Bellman-Zadeh, 
TOPSIS, and LINMAP decision making techniques for the Stirling 
cryogenic refrigerator and outcomes illustrated by Tyagi [9,48 . 

7. Error analysis 

To perform error analysis on the mentioned approaches, mean 
absolute percentage error (MAPE) was utilized properly. For this 
purpose, 30 different executions were utilized to gain optimum 
solutions by mentioned decision making techniques (that include 
TOPSIS, LINMAP and Bellman-Zadeh). The results obtained from 
error analysis are demonstrated in Table 2. According the MAPE 
values reported in able 2, LINMAP approach has a better solution 
for P and R L in contrast with other decision making methods; how¬ 
ever, fuzzy method has a great solution for COP draw in parallel 
with the rest decision making techniques. 

8. Conclusions 

In this study, to calculate input power, the cooling load ( R L ) and 
the coefficient of performance (COP) of the Stirling cryogenic 
refrigerator, FTT approach was employed. The aforementioned 
parameters of the Stirling cryogenic refrigerator were simulta¬ 
neously involved throughout multi-objective optimization process 
that includes following decision variables: 1 - The effectiveness of 
the cold-side heat exchanger (e L ), 2 - the hot-side heat exchanger’s 
effectiveness ( s H ), 3 - heat sink’s heat capacitance rate (C L ), 4 - heat 
source’s heat capacitance rate ( C H ), 5 - cold side’s temperature (T c ) 
and 6 - temperatures ratio (x). It is worth to mention that NSGA-II 


algorithm was also executed to evolve the performance of the mul¬ 
ti objective evolutionary algorithm used. After running mentioned 
optimization algorithm, the Pareto optimal frontier throughout 
objectives space was calculated. On the basis of the Pareto optimal 
frontier gained, an optimum solution was acquired by the three 
decision making approaches introduced previously that catego¬ 
rized in fuzzy, LINMAP and TOPSIS methods. Throughout this com¬ 
munication, it is proved this fact that the multi objective 
optimization method eventuates to favorable design of the Stirling 
cryogenic refrigerator draw a parallel with the single objective 
optimization method that conducted in previous literature [9,48 . 
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